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Abstract 

We construct a holographic dark energy model in a braneworld setup that gravity 
is induced on the brane embedded in a bulk with Gauss-Bonnet curvature term. We 
include possible modification of the induced gravity and its coupling with a canonical 
scalar field on the brane. Through a perturbational approach to calculate the effective 
gravitation constant on the brane, we examine the outcome of this model as a candidate 
for holographic dark energy. 
PACS: 04.50.-h, 98.80.-k, 95.36,+x 
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1 Introduction 



It is observationally confirmed from supernovae distance-redshift data, the microwave back- 
ground radiation, the large scale structure, weak lensing and baryon oscillations that the 
current expansion of the universe is accelerating [1]. There are several approaches to explain 
this late-time accelerated expansion of the universe. One of these approaches, is to introduce 
some sort of unknown energy component ( the dark energy) which has negative pressure (See 
[2] and references therein). The simplest candidate in this respect is a cosmological constant 
in the framework of the general relativity. However, huge amount of fine-tuning, no dy- 
namical behavior and also unknown origin of emergence make its unfavorable [3]. Another 
alternative is a dynamical dark energy: the cosmological constant puzzles may be better in- 
terpreted by assuming that the vacuum energy is canceled to exactly zero by some unknown 
mechanism and introducing a dark energy component with a dynamically variable equation 
of state [2]. Nevertheless, the main problem with the dark energy is that its nature and 
cosmological origin are still obscure at present. Other alternatives to accommodate present 
accelerated expansion of the universe are modified gravity [4] and some braneworld scenarios 
such as the Dvali-Gabadadze-Porrati (DGP) scenario [5]. Here, we are interested in to probe 
the nature of dark energy in the context of holographic dark energy models. The idea of holo- 
graphic dark energy comes from quantum gravity considerations [6]. This model proposes 
that if pa is the quantum zero-point energy density caused by a short distance cut-off, the 
total energy in a region of size L should not exceed the mass of a black hole of the same size. 
The holographic dark energy density corresponds to a dynamical cosmological constant and 
in this respect, the standard General Relativity should be modified by some gravitational 
terms which became relevant at present accelerating universe. Modified gravity provides the 
natural gravitational alternative for dark energy. Moreover, modified gravity presents a nat- 
ural unification of the early time inflation and late-time acceleration due to different role of 
gravitational terms relevant at small and at large curvature [4,7]. Also modified gravity may 
naturally describe the transition from non-phantom phase to phantom one without necessity 
to introduce the exotic matter. f(R) gravity is one successful attempt in this direction. 
Another theory proposed as gravitational dark energy is the scalar-Gauss-Bonnet gravity 
which is closely related with low-order string effective action [8]. The possibility to extend 
such consideration to third order (curvature cubic) terms in low-energy string effective action 
exists too [9]. Moreover, one can develop the reconstruction method for such theories [10]. 
It has been demonstrated that some scalar- Gauss-Bonnet gravities may be compatible with 
the known classical history of the universe expansion [7]. 

In this paper we propose a unified /(-R)-Gauss-Bonnet gravity with non-minimal coupling 
of the scalar field to f(R). We proceed a holographic dark energy approach to examine 
cosmological dynamics in this setup and we show that this model accounts for phantom-like 
behavior. This model presents a smooth crossing of the phantom divide line by its equation 
of state parameter and this crossing occurs in the same way as is supported observationally, 
that is, from above -1 to its below. 

2 The Setup 

The action of our modified GBIG model contains the Gauss-Bonnet term in the bulk and 
modified induced gravity term on the brane where a scalar field non-minimally coupled to 
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induced gravity is present on the brane 
1 
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where /?(> 0) is the GB coupling constant, k 2 is the five dimensional gravitational constant, 
A is the brane tension, A 5 is the bulk cosmological constant and K is the trace of the mean 
extrinsic curvature of the brane. 

To obtain cosmological dynamics in this setup, we use the following line-element 

ds 2 = q^dx^dx" + b 2 (y, t)dy 2 = -n 2 (y, t)dt 2 + a 2 (y, t)^ ij dx i dx j + b 2 (y, t)dy 2 , (2) 

where 7^ is a maximally symmetric 3-dimensional metric defined as 7^ = Sij + k and 
k = —1,0, 1 parameterizes the spatial curvature. The generalized cosmological dynamics of 
this setup is given by the following Friedmann equation ( see [11] for details) 
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where a dash denotes ^ and a dot marks J^. The energy-density corresponding to the 
non-minimally coupled scalar field is as follows 

1 
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where y is the coordinate of the fifth dimension and the brane is located at y — 0. To proceed 
further, we set n 2 (y, t)\ y= o = 1. Now we solve analytically the friedmann equation (3). It is 
convenient to introduce the dimensionless variables ( see for instance [12]) 
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where r = ^ is the DGP crossover scale and by definition 
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where 



Rf'(R) - f(R) - 6HRf"(R) 
The Friedmann equation (3) now takes the following compact form 

H 3 + H 2 + bH - p = 0. 



(9) 



The number of real roots of this equation is determined by the sign of the discriminant 
function N defined as [13] 

N = Q 3 + R 2 (10) 

where 
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If N > 0, then there is a unique real solution. If N < 0, there are 3 real solutions. Finally, 
if iV = 0, all roots are real and at least two are equal. In which follows, we consider the case 
with N > and in this case the unique real solution is given as follows [12,13] 
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Therefore, we achieve the following solution for the Friedmann equation (3) 
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(13) 



(14) 



(15) 



After obtaining a solution of the Friedmann equation, in the next section we try to 
provide the basis of a holographic dark energy in this setup. 



3 Effective Gravitational Constant 



Now we derive the gravitational constant within a perturbation theory to construct our 
holographic dark energy scenario. Firstly, we analyze the weak field limit of the model 
presented in the previous section within the slow time- variation approximation and at small 
scales with respect to the horizons size. It is convenient to use the longitudinal gauge metric 
written directly in terms of ( = log a 



ds 2 = e 2 < 



(i + 2*; 



a 2 H 2 



+ (1 - 2T)dx i dx i 



(16) 



The field * is the Newtonian potential and T is the leading-order, spatial post-Newtonian 
correction which will permit us to calculate the stress-anisotropy. For the longitudinal post- 
Newtonian limit to be satisfied, we require AT >> a 2 H 2 x (T, T, T) and similarly for other 
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gradient terms [14]. For a plane wave perturbation with wavelength A, we see that H 2 T is 
much smaller than AT when A << jj . The requirement that T is also negligible implies 
the condition 

^ll « J— (17) 
d( (A#) 2 ' 1 } 

which holds if the condition A << is satisfied for perturbation growth. This argument 

can be applied for T, T and Sip. Now using the metric (16), we can rewrite the and 
(ij) components of the gravitational field equations and the scalar field equation as follows 
[14] 
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where H is given by (15) and we have defined the following quantities 
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We also include the (Oj) component of the gravitation field equations in the same limit 

3 - a 2 H 2 9 n m = (B- A)AV + - A AT - B***. 

2 ip (p 

The energy-momentum conservation equations are [14] 

S = -9 

and 

H A^ 
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where 5 = Sp/p is the matter density contrast and 9 is the divergence of the matter peculiar 
velocity field. Finally, Poissons equation in this setup is as follows 
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(29) 
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This equation ( which its derivation can be obtained in [14] with more details) is the basis 
of our forthcoming arguments. After calculation of the effective gravitational constant, in 
the next section we construct a holographic dark energy model in the normal branch of 
this DGP-inspired modified GBIG scenario. We note that since the normal DGP branch 
is ghost-free and cosmologically stable, there is no ghost instability in our modified GBIG 
scenario too. 

4 Holographic Dark Energy in Modified GBIG Sce- 
nario 

4.1 General Formalism 

To study holographic dark energy model in the modified GBIG scenario, we first present a 
brief overview of the holographic dark energy model [6]. It is well-known that the mass of a 
spherical and uncharged D-dimensional black hole is related to its Schwarzschild radius by 



where the D-dimensional Planck mass, M D , is related to the D-dimensional gravitational 



constant G D and the usual 4-dimensional Planck mass through M D = G D D 2 and M 2 = 

M®~ 2 Vd-a with Vd-4 the volume of the extra-dimensional space. If pad is the bulk vacuum 
energy, then application of the holographic dark energy proposal in the bulk gives 



where V(S D 2 ) is the volume of the maximal hypersphere in a D-dimensional spacetime 
given by V(S D ~ 2 ) = Acr ' 1 . A d is defined as 



[15] 
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for D — 1 being even or odd respectively. By saturating inequality (27), introducing L as a 
suitable large distance ( IR cutoff) and c 2 as a numerical factor, the corresponding vacuum 
energy as a holographic dark energy is given by [15] 

Pad = c 2 (^M D ) D ~ Z M D A- D l D ( ~\ l- 2 . (33) 

8T {—) 

Using this expression, one can calculate the corresponding pressure via continuity equation 
and then the equation of state parameter of the holographic dark energy defined as wa = ^ 
can be obtained directly. 

Now we use this formalism in our GBIG setup. There are alternative possibilities to 
choose L ( see [15]). Here and in our forthcoming arguments we choose the IR cut-off, L , 
to be the crossover scale r which is related to the present Hubble radius via r = 1.2QHQ 1 ( 
see for instance the paper by Lue in Ref. [5] and [17]) where H is given by (15). In this case 
the effective holographic dark energy density is defined as follows 

*" = 8^ (34) 
where H and G e // are given by (15) and (30) respectively. So we find 

3H 2 f A 2 E + 2ABD - B 2 C^ 
P ^ = 8^G ( W^CE ) • (35) 

Using the conservation equation 

Peff + 3H(p eff + p eff ) = 0, (36) 

we can deduce 



Peff — H(D 2 + CE)~ 



(2DD + CE + CE) (A 2 E + 2ABD - B 2 C) - (D 2 + CE) (2AAE + A 2 E 
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Hence, the equation of state parameter of the model is given as follows 

1 (2DD + CE + CE)(A 2 E + 2ABD - B 2 C) - (D 2 + CE){2AAE + A 2 E + 2ABD + 2ABD + 2ABD - 2BBC - B 2 C) 2H + AH 



iff 



3H (A 2 E + 2ABD - B 2 C)(D 2 + CE) 3H 2 

(38) 

Now to proceed further, we should specify the form of f(R). In which follows we present a 
relatively general example by adopting the Hu-Sawicki model. 



4.2 An Example 

As an example, we consider the following observationally suitable form of f(R) ( the Hu- 
Sawicki model [16]) 

« s »= s -^ < 39 > 
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where for an spatially flat FRW type universe, 

R=-Q(H + 2H 2 ) (40) 

and both and R c are free positive parameters. We also adopt the ansatz with a = a^f and 
if = <^o£~ M (with v and /i positive constants) and a((p) — |(1 — £y? 2 ) with £ a positive non- 
minimal coupling. We set V(ip) = Vo<^ 2 and we translate all of our cosmological dynamics 
equations in terms of the redshift using the relation between a and z as 1 + z = 211 . One 
can use these ansatz in the Friedmann, Klein-Gordon and the conservation equations to find 
constraints on the values of v and \i ( and also £) to have an accelerating phase of expansion 
at late-time ( to see a typical analysis in this direction see [17]). Based on such an analysis, 
in which follows we set v = 1.2 and fi = 0.5. The value adopted for v gives a late-time 
accelerating universe. Then we perform some numerical analysis of the model parameter 
space. We set £ = 0.15 which is close to the conformal coupling. The reason for adopting 
such a value of £ is motivated from constraint on £ from recent observations [18]. For our 
numerical purposes we have set Vq = 1 and (3 = 0.33. Figure 1 shows the behavior of the 
effective energy density, p e ff, versus the redshift. We see that the effective energy density 
grows with time and always p e ff > 0. This is a typical behavior of phantom-like dark energy 
density. Figure 2 shows the behavior of 1 + u e ff versus the redshift. With the choice of 
the model parameters as we have adopted here, the universe enters the phantom phase at 
z ~ 0.25 close to the observationally supported value. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



Figure 1: Variation of the effective dark energy density versus the redshift. 
The deceleration parameter defined as 
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Figure 2: Variation of 1 + uj e ff versus the redshift. 



in our setup takes the following form 
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Figure 3 shows the behavior of q versus z. In this model, the universe has entered an 
accelerating phase in the past at z ~ 0.84. 



q 




Figure 3: Variation of the deceleration parameter versus the redshift. 



Variation of H with cosmic time or redshift gives another part of important information 
about the cosmology of this model. We deduce the relation for variation of H versus the 
cosmic time as follows 
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Figure 4 shows the variation of versus the redshift. Since H < always, the model 
universe described here will not acquire super-acceleration and big-rip singularity in the 
future. 



\dot{ff} 
Hi 




Figure 4: Variation of jL versus the redshift. 
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This example shows how a modified GBIG scenario has the potential to give a phantom- 
like behavior in a holographic viewpoint. This behavior is realized without introducing any 
phantom field in this setup. In fact a combination of the curvature effect and the non-minimal 
coupling provides the suitable framework for realization of this phantom mimicry. 



5 Summary 

The model presented here, contains the Gauss-Bonnet term as the UV sector of the theory, 
while the Induced Gravity effect completes the IR side of the model. The induced gravity 
on the brane is modified in the spirit of f(R) gravity which itself provides the facility to self- 
accelerate even the normal, ghost-free branch of the DGP-inspired model [11,19]. The model 
also considers a non-minimal coupling between the canonical scalar field on the brane and 
the modified induced gravity. This is a general framework for treating dark energy problem 
and other alternative scenarios can be regarded as subclasses of this general model. We have 
studied the cosmological dynamics in this generalized braneworld setup within a holographic 
point of view. By adopting a relatively general ansatz for f(R) gravity (Hu-Sawicki model) 
on the brane, we have shown that the modified GBIG scenario presented here realizes the 
phantom-like behavior: the effective energy density increases with cosmic time and the ef- 
fective equation of state parameter crosses the phantom divide line smoothly in the same 
way as observations suggest: from quintessence to the phantom phase. We have studied 
also the holographic nature of the dark energy in this setup by calculation of the effective 
gravitational constant via a perturbational approach. In this model, H < always, and 
therefore the model universe described here will not acquire super-acceleration and big-rip 
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singularity in the future. 
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